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Abstract 

The Lagrangian for the motion of n well-separated BPS monopoles is calculated, by 
treating the monopoles as point particles with magnetic, electric and scalar charges. It 
can be reinterpreted as the Lagrangian for geodesic motion on the asymptotic region 
of the n-monopole moduli space, thereby determining the asymptotic metric on the 
moduli space. The metric is hyperkahler, and is an explicit example of a type of metric 
considered previously. 
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The moduli space for n Bogomolny-Prasad-Sommerfield (BPS) monopoles, M n , is An- 
dimensional, and has a natural hyperkahler metric [1]. It has the form M n = M 3 x 
(S 1 x M°)/Z n , where M° is simply connected and is the moduli space of "strongly 
centred" monopoles [2]. The n-fold cover of M n , M 3 x S 1 x M°, is metrically a direct 
product, with M 3 x S 1 being flat. The geodesies on M n accurately model low energy 
monopole dynamics [3,4]. 

It is generally understood that when n monopoles are well-separated in space (R 3 ), the 
4n coordinates on M n can be thought of as positions and phases of n indistinguishable 
particles. Bielawski has shown that to a first approximation, and working locally, 
the metric is flat in these coordinates, and he has identified the coordinates with the 
parameters occuring in the rational maps which are associated with monopoles [5]. 
However, the detailed topology of this outer region of moduli space has not been clearly 
elucidated, neither is the metric known accurately. It is known that if n > 1, it is not 
just the product of a configuration space of n point particles by an n-torus. In this 
paper we shall calculate the metric on this outer region of moduli space, and make 
some remarks on the topology. 

In the simplest case, n = 1, the moduli space M\ = R 3 x S 1 and the metric is flat. The 
geodesic motion is at constant velocity in M 3 , and constant angular velocity 6 on the 5" 1 - 
factor. Physically this describes a monopole with an electric charge (a dyon) in uniform 
motion; its magnetic charge g is fixed and its electric charge is proportional to 9. The 
geodesic motion on M n , according to Bielawski's metric, describes the independent 
uniform motion of n dyons in M 3 , each with magnetic charge g and a constant electric 
charge. 

Now, the true motion of monopoles, even if they are well-separated, is more interesting. 
The 2-monopole moduli space M 2 is known in detail, and its metric has been determined 
by Atiyah and Hitchin [6]. The asymptotic metric on M® , when the two monopoles are 
well-separated, is obtained by neglecting all terms in the Atiyah-Hitchin metric which 
are exponentially small, and it has a simple algebraic form. It is the self-dual euclidean 
Taub-NUT metric [7] with a negative mass parameter. Geodesies on Taub-NUT space 
describe two monopoles or two dyons interacting via Coulomb-like forces. 

The Taub-NUT metric has a U(l) symmetry not possessed by the Atiyah-Hitchin 
metric. This implies that well-separated dyons have individually conserved electric 
charges. In contrast, the general geodesic motion on Mi allows for electric charge 
exchange (as well as momentum exchange) in a close encounter of two monopoles or 
dyons. 
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Shortly after Atiyah and Hitchin obtained their metric, one of the present authors 
showed that its asymptotic form can be obtained from a physical calculation of the 
dynamics of two well-separated dyons [8]. It suffices to consider the dyons as point 
particles, each with a magnetic charge, electric charge and scalar charge. The equations 
of motion for the two dyons, assuming that their speeds are modest and that their 
electric charges are much less than their magnetic charges, are found to be equivalent 
to the equations for geodesic motion on Taub-NUT space, together with a centre of 
mass motion. The scalar charge can be understood from the exact solution of the field 
equations for a single dyon. For a magnetic charge g and electric charge q, the Higgs 
field far from the dyon centre has magnitude 



v - {9 +Q )2 + 0(e- 8 ™ r / 5 ), (1) 



4nr 



where v is the vacuum expectation value of the Higgs field. This may be interpreted 
as saying that the dyon has scalar charge (g 2 + q 2 )^ . 

In this paper we shall repeat this calculation, but for n dyons. Our strategy is to find 
the Lagrangian describing the dyons' motion in IR 3 , with the (constant) electric charges 
as parameters, and then to reinterpret the n electric charges as arising from motion 
on n circles associated with the n monopoles. In this way, the Lagrangian is shown 
to be equivalent to the Lagrangian for geodesic motion on a particular n-torus bundle 
over the configuration space of n monopole positions, with respect to a metric whose 
explicit form we give. 

The metric we obtain has some beautiful features, and its explicit form is quite simple. 
It admits an isometric n-torus action, as well as the action of the Euclidean group 
on M 3 . The equations of motion for the dyons, to which it gives rise, are Galilean 
invariant, which is not necessarily the case for a general moduli space Lagrangian. 
Finally it is hyperkahler. Hyperkahler metrics on 4n-dimensional spaces have been 
studied previously. A number of authors have presented formulae for such metrics 
involving matrices subject to some differential constraints. Our metric is in the class 
previously considered, and one may easily show that it satisfies the required constraints. 



Consider n dyons, all with the same magnetic charge g. Let the ith dyon have electric 

i 2 
charge q t and scalar charge (g 2 + q 2 ) 2 , which we shall later approximate by g+ its 



rest mass is rrii = v(g 2 + q 2 ) 2 . Suppose the dyons have positions x^ and velocities Vj. 
Denote their separations Xj — Xj by rji, and set rji = \rji\. The Lagrangian for the 
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motion of the nth dyon in the background of dyons 1, . . . , n — 1 is 

Ln = (-m n + (g 2 + ql)$<ty (1 - v*)* 

+ q n v n • A - g n A ( 2 ) 
+ #v n • A-5rA . 

Here is the scalar field due to dyons 1, . . . , n — 1; its effect is to modify the rest mass 
of dyon n, with the coefficient (g 2 + qfy^ being the scalar charge of dyon n. A, Aq 
are the vector and scalar Maxwell potentials due to all but the nth dyon, and they 
couple to the electric charge of the nth dyon, q n . The magnetic charge of the nth dyon, 
g, couples to the dual Maxwell vector and scalar potentials A, Aq produced by dyons 
1, . . . , n — 1. These potentials are defined so that 

V x A = B 
-VA - A = E 

(3) 

V x A = -E 
-Vl - A = B, 

where E and B are the electric and magnetic fields at the nth dyon, due to the other 
dyons. 

The fields <f>, A, Aq, A and Aq are linear combinations of the fields produced by the n — 1 
contributing dyons. The contribution from dyon 1 is given in terms of the Lienard- 
Wiechert potentials at x n due to the moving particle at xi. The scalar field due to 
dyon 1 at x n is 

*=^±^(l-vS». (4) 

Eqn. (4) is the Lorentz scalar version of a Lienard-Wiechert potential, where s n \ = 

( r ni — ( r m x v i) 2 + °( v i)) 2 • Now the leading terms in the Lagrangian will turn out 
to be of order v 2 /r, and this is the order in v to which we wish to work, so it is 
sufficient to approximate s n i by r n \. Expanding out the square roots, and keeping 
terms quadratic in velocities and electric charges, we find 

,=_£_fi + i?-4y ( 5) 



Anr nl V 2g 2 

To write down the Maxwell and dual Maxwell potentials, we introduce a Dirac vector 
potential w(y). This is defined to satisfy 

Vxw = 4 , (6) 



and w(y) = w(— y), but we do not need to specify a particular gauge. We denote by 
Wji the potential w(xj — Xj) = w(r ji); it has the symmetry property Wji = Wy and 
satisfies Vj x = V^l/r^). The potentials A, A 0: A and A produced by dyon 1, 
and evaluated at x n are 



= ~ fw nl • Vl (8) 

io = j^— + | L w nl -v 1 . (10) 
47rr n i 4/T 

These are boosted versions of the familiar Coulomb and Dirac potentials of a dyon at 
rest, and we have again approximated s n \ by r n \. Substituting (5) and (7) - (10) in 
(2), and keeping terms of order v 2 , gv and q 2 , we obtain 

L n = -m n + -m n v 2 - (v n - vi) 2 

nl (11) 

- j-(Qn ~ 9i)(v n - vi) • w n i + — *— (q n - q x ) 2 . 
4tt 87rr n i 

Eqn. (11) includes the contribution from dyon 1 to the Lagrangian for dyon n. Adding 
the contributions of all dyons 1, . . .n — 1 gives 

L n = -mv n - — ^ — — ^(q n - 9! )(v n - Vi) ■ w m 

(12) 

| 1 ft) 2 

(We have dropped the constant — m n , which has no effect on the equation of motion, 
and in the kinetic term replaced m n by m = vg, which is accurate enough to the order 
we are working.) It can be seen that the interaction terms are quite symmetric between 
dyons i and n. It follows, that if we symmetrize (12) between all the dyons, we will 
obtain a Lagrangian whose Euler-Lagrange equations are the equations of motion for 
all the dyons. The symmetrized Lagrangian is 

r 1 2 9 2 ( v j- v *) 2 9 sr^ i w x 

1 = 1^1™,-- ^ ^~ 47r 2. (^"ftKvi-vO-w^ 

i=l l<i<j<n J l<i<j<n 

- - " " (13) 

(ft - ft) 2 



87r r,i 
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(Here we have used the symmetry property Wji = Wy.) The kinetic term may be 
rewritten as 

n 1 1 1 

i=l l<i<j<n 

which shows that the sum of the dyon velocities completely decouples from the relative 
velocities. It follows that ^(vi + . . . + v n ) is a conserved quantity; it can be identified 
as the centre of mass velocity. Furthermore, the relative motion is unaffected by a 
Galilean transformation Vj — > Vj + V, for any fixed V. There is also invariance under 
the Euclidean group of translations and rotations. 

The equation of motion for dyon j is 

9 2 \- Vj - Vi g 2 ^ 2( Vj - v^Tji ■ (\j - v,) - (vj - Vi) 2 ^ 



mv 



4tt ^-f r^-i ' 8tt ^ 



(15) 



The last term in (15), an electric Coulomb force, depends on the square of the difference 
of electric charges rather than on their products. This is one of the consequences of 
the scalar interaction. Note that there is a static solution of the equations of motion 
(vi = : 1 < i < n) with all electric charges equal. This is to be expected, as there 
are solutions of the Bogomolny equations representing n monopoles at rest, and these 
can be transformed in a simple way (discovered by Julia and Zee [9]) into a stationary 
solution representing n dyons each with the same electric charge. 

The Lagrangian (13) is defined on the 3n-dimensional configuration space (7 n (M 3 ) = 
(R 3 ) n — A with coordinates {x^ : 1 < i < n}, where A represents the subspace of 
(R 3 ) n where two or more monopole positions coincide. The electric charges q t are 
constant parameters. The Lagrangian is not purely kinetic, because of the terms linear 
in velocity (the electric-magnetic couplings), and because of the Coulomb terms. But 
if we could interpret the electric charges as velocities, as in Kaluza-Klein theory, then 
the whole Lagrangian would be quadratic in velocities. We shall now show that this is 
indeed possible. 

Consider a 4n-dimensional manifold E n , which is a T n (n-torus) bundle over C n (M. 3 ) 
with local coordinates {xj,#i}. 9i is a phase angle associated with the ith monopole. 
Endow this bundle E n with a T n -invariant metric, so that the purely kinetic Lagrangian 
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for motion on E n is 



£ = \9ijVi ■ Vj + TjhijiOi + W ik ■ v k )(9j + Wji ■ v/), (16) 

where gij and hij are symmetric matrices. (From now on we shall be using the sum- 
mation convention unless we say to the contrary.) Invariance under the torus action 
requires that the matrices gij,hij and Wij depend only on the 3n coordinates {x^}. 
The Euler-Lagrange equations for C admit the n constants of motion 

q i = Kh ij (e j + W jl -Vi). (17) 

For a suitable value of the constant qi may be identified with the electric charge of 
monopole i. 

It would not be correct to substitute for these constants of motion in the second term 
of £, and then calculate the remaining equations of motion. However, it is correct to 
reduce the Lagrangian C to 

1 1 1 

where /i lJ is the inverse of hij. The equations of motion given by C and £ e ff are the 
same. It follows that we can determine gij,hij and by requiring that £ e ff and 
the Lagrangian L, eqn. (13), give the same equations of motion. This almost means 
the same thing as identifying £ e ff and L. The matrix gij must be chosen to have 
components 

9jj = m ~ ~a~ — ( no sum over ^ 

9 « = f^- (20) 



and Wjj must have components 



Wjj = — ^ ^ Wij (no sum over j) (21) 

Wy = f£w y (i^j). (22) 

The symmetry properties of the Dirac potentials imply that Wjj is a symmetric matrix. 
Simply identifying £ e ff and L would give a matrix /i* J with no inverse. But we may 
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add any constant matrix to /i lJ without changing the equations of motion coming from 

. . 2 

jCeff- Taking advantage of this we see that a satisfactory choice is h lJ = ^jgij- 

The value of k can now be fixed. The expressions (21) and (22) for Wy could lead to 
singularities in the Lagrangian (16), because the Dirac potentials have the usual Dirac 
string singularities. There are two Dirac strings in the potential w(y), because of the 
symmetry of the potential, and each carries "flux" 2n. These strings can be gauged 
away if k = 4ir/g, and if the angles 9i have the usual range [0, 2ir]. 

To obtain the Lagrangian C in its simplest mathematical form, it is convenient to 
choose units so that g = 4-7r,m = 4-zr (these are in any case true for BPS monopoles 
in the natural units [1]) and to remove an overall factor of 4tv. C is the purely kinetic 
Lagrangian for motion on E n with metric 

ds 2 = gijditi ■ ditj + g~^(dQi + W ik ■ dx k )(d9j + Wji ■ dx t ) (23) 

where, now, 

gjj = 1 — ^ — (no sum over j) (24) 





1 


9zj 




Tij 




= -£ 






Wy 


= Wy 



(i^j) (25) 
Wy (no sum over j) (26) 

(i*3)- (27) 

Notice that if all terms which vary inversely with monopole separations are neglected, 
then the metric (23) reduces to Bielawski's flat metric. 

It can be easily verified that the matrices <7y and Wy satisfy 



^ ~ r b <^_wa _ ^abc & 

«5 



^-M w& = e ° ie ^' (28) 



d 

= Q^adik, (29) 

where xf and denote the components of Xj and Wy. These conditions, (28) & 
(29), were shown by Pedersen and Poon [10], and Papadopoulos and Townsend [11], 
following earlier work by Hitchin et al. [12] to be the necessary conditions for the metric 

(23) to be hyperkahler. These authors gave no explicit solution. Our explicit solution 

(24) - (27) is surprisingly rather simple and symmetric. In fact, eqns. (28) & (29) may 
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be solved locally by a rather more general ansatz than (24) - (27). It suffices to replace 
l/rij by a function iJy(|xj — Xj|) which is harmonic in both arguments. The Dirac 
potential Wy is then chosen to satisfy Vi x Wy = ViHij (no sum over i). 

To obtain the asymptotic region of the n-monopole moduli space M n , one must quotient 
E n by the permutation group S n , which acts by permuting the positions and phases 
of the monopoles. This is because the monopoles are unordered, or indistinguishable. 
One also requires ^> 1, for all i ^ j. Thus in the case of two monopoles, where the 
metric calculated using the method above is the Taub-NUT metric times a flat factor, 
one must identify under the interchange of (xi,#i) and (x2,#2) 5 and have T\i ^> 1, in 
order to reproduce the asymptotic form of the moduli space. 

The three Kahler forms on E n are 



It follows that the torus action is triholomorphic 



and thus that the generators of the torus action, i.e. the Killing fields are Hamil- 
tonian vector fields with respect to all three symplectic forms u a . The three moment 
maps associated with -J=^- are the three coordinate functions xf . By contrast the rota- 
tion group SO (3) does not act triholomorphically; the forms u a transform as a triplet 
under its action. 

Using the moment maps xf (all of which Poisson commute with respect to a fixed 
symplectic structure) one may obtain lower-dimensional hyperkahler metrics by taking 
the hyperkahler quotient. This amounts to setting k of the x^'s to constants and then 
projecting the metric orthogonally to the k vector fields In this way one obtains a 
new 4(n — /c)-dimensional hyperkahler metric. One may use this freedom to freeze all 
but one of the coordinates Xj. The result is a four-dimensional multi-centre metric of the 
form discussed in [13]. If one wishes, one may fix not the x^'s but linear combinations 
of them. A case of particular interest is when one reduces using the position of the 
centre of mass 



This amounts to taking out the centre of mass motion. In fact, because of the special 
form of the metric (i.e. eqns. (24) - (27)), it is easy to see directly that it splits as a 




(30) 




(32) 
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metric product of the centre of mass manifold, with coordinates X and 9 = $i, and 
a moduli space of centred monopoles. 

Let us now turn to the topology of E n . The space C n (R 3 ) is the space of n distinct 
ordered points in R 3 . It is simply connected and its second homology group H2(C n (M. 3 )) 
is |n(n — 1) -dimensional and torsion free. A homology basis is provided by the ^n(n — 1) 
2-spheres Sfj : |xj — Xj| = = constant, = constant, j ^ k ^ i. To specify a 

T n bundle over C n (R 3 ) it suffices to specify its behaviour over these 2-spheres. A T n - 
invariant metric gives rise to a connection on the bundle E n if we define the horizontal 
subspaces of the tangent space of E n to be orthogonal to the fibre directions with 
respect to that metric. In our case this means that the connection one-form associated 
to the ith generator of T n is given by d9i + • <ix/j . Evidently, restricted to Sfj , this 
reduces to the standard Dirac monopole connection over 5" 2 with unit magnetic charge. 
This specifies the T n bundle on (7 n (R 3 ). 

In conclusion we have shown that the equations of motion for n well-separated BPS 
monopoles or dyons can be obtained from the purely kinetic Lagrangian for motion on 
a 4n-dimensional manifold E n /S n , where E n is an n-torus bundle whose metric (23) is 
hyperkahler, and S n is the permutation group. The particle motion in R 3 corresponds 
to geodesic motion on E n /S n . The space E n / S n is fibred over C n (R 3 ), the configuration 
space of n indistinguishable particles in R 3 . The particles must be at distinct locations, 
as there are singularities if they coincide. E n /S n is therefore incomplete. Moreover 
if the particles are sufficiently close our metric is no longer positive definite. This 
happens because the "mass" parameters, i.e. the coefficients of the second term in 
(24), are negative. 

On the other hand it is known that the low energy dynamics of n BPS monopoles 
or dyons at arbitrary separation is accurately modelled by geodesic motion on the 
moduli space M n of static solutions to the Bogomolny equations, and M n is complete 
and thus has no singularities. Thus E n /S n only gives the asymptotic form of M n , 
for well-separated monopoles. In the case n = 2, it is known that E n /S n and M n 
differ by an exponentially small amount as the separation gets large. It is natural 
therefore to conjecture that E n /S n and M n are exponentially close, for any n, if all 
pairs of monopoles are well-separated. This would imply that in the scattering of dyons, 
there would only be exponentially small electric charge exchange, provided the dyons 
remained well-separated, although momentum exchange falls off with a power of the 
separation. 

Mathematically one may consider a version of the metric on E n for which the mass 
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parameters are positive. This metric is everywhere positive definite. For suitably 
chosen values of the mass parameters it is also complete, as the subspace where pairs 
of points in M 3 coincide corresponds to the fixed point set of one of the generators of 
the torus action. The simplest case is n = 2 which gives the usual Taub-NUT metric 
with positive mass. We believe that the global structure of these positive definite 
metrics would repay further study in view of their application as target spaces for 
supersymmetric sigma-models. 

We have made no progress with solving the equations of motion for n well-separated 
dyons. In the n = 2 case, we found enough constants of motion to determine the asymp- 
totic 2-dyon scattering and bound orbits explicitly [14]. It would be very interesting 
if the n-body problem was tractable. We have however found a simple expression for 
the total angular momentum which includes the contribution from the interaction of 
electric and magnetic charges, and which generalizes that found in [14]. It is 



J = 47r $^#u x * x v j - - — , (33) 



Til 

1,3 i<j 



or equivalent ly, 



J = 4^x ! xv t - £ 4 7 r(^-x J )x(v,-v J ) + (x,-x J )fa- 9j )_ ^ 



i=l l<i<j<n ^ 



It is interesting to study the quantized dynamics of monopoles, and especially the 
dynamics of monopoles in super- Yang-Mills theory. Sen has recently argued that there 
are certain quantum bound states of n monopoles with zero energy in such a theory 
[15]. To verify Sen's conjectures one first has to understand the cohomology of the 
n-monopole moduli space M n , which depends largely on the properties of M n when 
all monopoles are close together. However to check in detail whether the cohomology 
classes are directly related to normalizable quantum states requires a careful analysis of 
the metric on M n for well-separated monopoles, and here our results may be useful. In 
the simplest case of two monopoles, the bound state wavefunction decays exponentially 
with separation and it may be that for n monopoles, the wavefunction similarly decays 
if any one monopole becomes separated from the rest. 

Finally we remark that if we set the electric charges to zero in (13) we obtain a metric 
on (R 3 ) n which is up to rescaling almost identical to that obtained by Shiraishi [16] for 
n stringy black holes. The difference is that in Shiraishi's case the minus sign between 
the first two terms in (13) is replaced by a plus sign. This would seem to imply that 
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by adjoining extra angles to the moduli space of n stringy black holes one may obtain 
a hyperkahler moduli space. This may lead to a resolution of the puzzle discussed in 
[17]. 
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